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Abstract 

We analyse the radiative decay r/ — »• 7r"'"7r~7r°7 in the low-energy expansion of the Stan- 
dard Model. We employ the notion of "generalized bremsstrahlung" to take full advan- 
tage of the theoretical and experimental information on the corresponding non-radiative 
— > Stt decay. The direct emission amplitude of 0{p^) is due to one-loop diagrams with 
intermediate pions (isospin violating) and kaons (isospin conserving). The isospin con- 
serving contributions to direct emission, including vector meson exchange appearing at 
C(p^), are suppressed. 
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1. The decays — > Svr are forbidden in the limit of isospin conservation. Neglecting 
the small electromagnetic corrections [ |T]], the amplitudes are proportional to the isospin 
breaking mass difference m„ — m^. The leading-order amplitude in the low-energy ex- 
pansion of C(p^) [ 0] is known to receive large higher-order corrections, both at O(p^) [ 
§ and beyond [|,|5|]. 

The radiative decay i] 7r+7r^7r°7 is in principle an interesting channel. At lowest 
order p'^, the amplitude is pure bremsstrahlung. At next-to-leading order an additional 
contribution appears (direct emission) that is non-vanishing even in the isospin limit. 
Therefore, the direct emission amplitude carries in principle new information that is not 
accessible in ^ Stt decays. The notion of a direct emission amplitude is not unique 
except that it starts at 0{k) where k is the photon momentum. For instance, the so-called 
quadratic slope parameters of the non-radiative amplitude arising at 0{p'^) also generate 
a radiative amplitude of 0{k) that one may combine with the bremsstrahlung amplitude 
because it is also completely fixed by the non-radiative process. We have recently shown 
[ P that one can define a generalized bremsstrahlung (GB) amplitude for a generic radia- 
tive four- meson process that includes the effects of all local terms of 0{p'^) contributing 
to the non-radiative transition. 

The main advantages of the GB amplitude are: 

• Since all local contributions to the non-radiative amplitude of 0{p^) are included, 
the uncertainties in the corresponding low-energy constants do not propagate to the 
direct emission amplitude (defined here as the difference between the total and the 
GB amplitudes). 

• If there are substantial higher-order contributions beyond 0{p'^) in the non-radiative 
amplitude they can be included in the GB amplitude by using the experimentally 
measured non-radiative amplitude. For rf 7r"'"7r~7r°7, this is especially welcome 
because the unitarity corrections [ ^ |]] modify both rate and slope parameters of 
?7 ^ Stt substantially. Using the experimental values in the GB amplitude allows 
for a much more accurate determination of the total amplitude. 

The purpose of this letter is to calculate both GB and direct emission amplitudes for 
?7 7r+7r^7r°7 along the same lines as for K 37r7 [0]. We comment on the differences 
between the GB and the usual bremsstrahlung amplitudes and we discuss the relative 
importance of the main contributions to direct emission: pion loops (isospin violating), 
kaon loops and vector meson exchange (both isospin conserving). The suppression of the 
isospin conserving component of direct emission is explained. 
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The channel under consideration has already been studied in the framework of chiral 
perturbation theory by Bramon et al. [ |^] where also references to the earlier literature 
can be found. We will discuss the differences to Ref. [ |Sp as we go along. 

The upper limit B{r] 7r+7r^7r°7) < 6 x 10""^ quoted by the Particle Data Group [^] 
refers to direct emission only [ |T^. The experimental situation will improve considerably 
in the near future. For instance, the KLOE experiment at the Frascati ^-factory [ [TT|] 
should collect more than 10^ r] per year. 



2. To evaluate the bremsstrahlung contribution to r] ^ 7r+7r^7r°7 we need to know the 
amplitude for ?7(]9^) — vr+(p+)7r^(]9_)7r°(po) • Neglecting electromagnetic corrections [ 
the amplitude can be written in the form [ ^ 



where 5 is a parameter of the lowest-order chiral Lagrangian [ |T2l] related to the quark 
condensate and F.„ = 92.4 MeV is the pion decay constant. The kinematical variables 
s, s±, So are defined as 

s = {Pr,- Pof , s± = - P±f , So = ^(s + s+ + s_) . (2) 

The function s±) vanishes to lowest order p^. At 0{p^) it receives both loop and 
counterterm contributions [ ^ . Higher-order effects due to vrvr rescattering are important 
and have been included in 5(s, s±) by way of dispersion relations [ 0, |]]. These higher- 
order corrections increase the rate of 0{p'^) by some 25 ^ 30 % and must be included for 
a reliable estimate of the bremsstrahlung amplitude. 

Experimental results are conventionally expressed in terms of the Dalitz variables x, y 
defined as 

2M^Q ' ^ 2M^Q ^ ^ ' ^ ' ^ ^ 

Q = M^- 2M^+ - M^o . 

Up to a normalization constant, the experimental Dalitz plot distribution is fitted by a 
function of the form [ ^ |T3]] 

A{x, yf = A{0, 0)2(1 + ay + by^ + cx^) (4) 

where A{x,y) corresponds to the decay amplitude ([T]). Charge conjugation invariance 
forbids a term linear in x. 
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a 


b 


c 


Experiment [Ell^-^ 


-1.22 ±0.07 


0.22 ±0.11 




Gasser and Leutwyler C(p^) ^ 


-1.33 


0.42 


0.08 


Kambor et al. (solution a) 


-1.16 


0.24 


0.09 


Kambor et al. (solution b) 


-1.16 


0.26 


0.10 



Table 1 : Experimental and theoretical values of the linear and quadratic slopes of t] ^ 
7r+7r^7r° defined in Eq. (Q). 

The present experimental and theoretical status of the parameters in is summarized 
in Table [l|. We do not need a value for A{0, 0) since we always normalize our results to 
the non-radiative decay. In this way, errors are substantially reduced. From Table |T], the 
importance of higher-order corrections is evident also for the slope parameters. For the 
numerical calculation, we will use the experimental values of a,b. Experiments have not 
been sensitive enough to extract the parameter c which is however relatively stable with 
respect to chiral corrections (we will take c = 0.10 for the numerics). 

The kinematics of the decay ?7(p,,) 7r+(p+)7r^(p_)7r°(po)7(fc) is specified by 
adding the variables 

ti = k-pi (i = r],+,-,0) (5) 

with 

t,j = t+ ± t_ ± to • 

Any three of the tj together with x and ?/ in form a set of independent variables. 
With CP conserved, there is only an electric transition amplitude that we write as 

A{r] 7r+7r-7r°7) = ee''{k)E^ (6) 

with 

fc^E^ = . 



Low's theorem [ IT^ relates the radiative amplitude to the corresponding non-radiative 



amplitude and their first derivatives with respect to the Dalitz variables up to 0{k). For 

a general 

dA{s,t) 



dA(s t) 

a general four-body amplitude A(s,t) with Mandelstam variables s, t, both — — and 

OS 



contribute to the Low amplitude. Since there are two neutral particles in our pro- 

ot 

cess we can choose variables and assign particle labels such that only one of the deriva- 
tives enters. With pi = — p^, p2 = Po, Ps = P- and p^ = p+ in the notation of Ref. [ §], 
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Low's theorem reads 



where 



X 



y 




(7) 



+ oik) 



(8) 



from now on. 

To lowest order p^, the radiative amplitude is completely given by the Low amplitude 
(^. In fact, since there is no x-dependence in the 77 ^ Svr amplitude of 0{p^) in only 
the non-derivative part in contributes. Starting at 0{p^), an x-dependence is generated 
that produces the quadratic slope term cx^ in (Q). 

However, one can do better than that. In order to account for all the local parts of 
0{p'^) in the non-radiative amplitude that contribute also to the radiative amplitude, a so- 
called generalized bremsstrahlung amplitude can be introduced [ ^. One major advantage 
of using the GB amplitude is that the remaining direct emission amplitude E'^ — Eq^ 
can only receive contributions from local terms of C(p^) that do not contribute to the 
non-radiative amplitude. For r] — > 7r+7r^7r°7, only the low-energy constant Lg [ [1^] 
could therefore appear in the direct emission amplitude. However, the corresponding 
counterterm does not contribute to — 7r+7r^7r°7 even for m„ 7^ m^. Thus, the one-loop 
contribution to direct emission is necessarily finite. 

The general formula for the GB amplitude of Ref. [ 

to 



simplifies in the present case 



-^GB 



+ 



73 
3 



{tr,P^ - top';, 



Po+Pn 



dA{x,y) 
dx 

{to + tr,) 



pt 



d'^A{x, y) 
dxdy 



{t-P+ - t+P^- 



d'^A{x, y) 
dx'^ 



+ 0{k) . 



If one uses the experimental amplitude as given by the Dalitz plot distribution 



(9) 
the 



last term in (^) involving 



d^A{x, y) 
dxdy 



will in fact not contribute. As already announced. 
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(MeV) 


rGB(^ — ^ TT + TT 7r*^7)/r(?7 ^ TT + TT 7r°) 


10-30 


(2.30 ±0.04) X 10-3 


30-50 


(5.99 ±0.10) X 10-^ 


50-70 


(1.85 ±0.04) X 10-^ 


70-90 


(4.47 ±0.11) X 10"^ 


> 90 


(5.00 ±0.14) X 10-6 



Table 2: Rates for V{ri ^ tt+tt 7r°7) with the GB amplitude (Q) for different bins in the 
photon energy E^, normalized to r(77 ^ n^n^-n^). 



we use for the slope parameters the experimental values [ ^ [T3|] a = —1.22 ± 0.07, 
b = 0.22 ±0.11 and the theoretical prediction [0] c = 0.10. The results for the rate 
normalized to T{t] ^ n^n^n^) are given in Table |^ for five bins in the photon energy E^ 
(in the r] rest frame). 

The relative branching ratios for E^ > 10 and 50 MeV, respectively are 



B{7] ^ n+n-n\,E^>10MeV)GB = (3.14 ± 0.05) x 10-^5(77 ^ tt+tt-ttO) 
B{7] ^ TT+TT-n\Ey>50MeV)GB = (2.35 ± 0.05) x 10-^5(7] ^ tt+tt-ttO) . 

(10) 

The errors given in both ( [T0| ) and Table || are due to the experimental errors of the slope 
parameters a, b. These errors would of course be much larger if we would not normalize 

to rij] 71^71^71^). 

We can now make a first comparison with the work of Ref. [ ||]. Bramon et al. con- 
structed a simple approximation to the Low amplitude (^. They dropped the derivative 
term in (0) and took instead the amplitude A(x, y) of 0{p^) [ ^ at the center of the Dalitz 
plot. In fact, they did not exactly use the amplitude of Ref. [ ^ but increased the coun- 
terterm amplitude to account for the discrepancy between the experimental rate and the 
predicted rate of 0{p^). With these assumptions, they obtain [ ISO 



TT + TT 7r°) 



Bij] ^ TT+TT-TT^; E^ > 10 MeV)bremsstrahlung = 2.81 X l{)-''B{'n 

B{r] ^ n+n-n\, E^ > 50 MeV)brem..trahiung = 1.85 x 10-*B{r] ^ n+n-7v^) 

(11) 

In spite of the rather drastic approximations made, this prediction is quite close to our 
result (ITD|) that is based on the GB amplitude (1^) and on experimental input for the slope 



parameters. Of course, it is difficult to assign an error to the prediction of Bramon et 
al. With the errors given in due to the experimental errors of the slope parameters, 
our prediction for B(r] n^n^n^'j; E^ > 10 MeV) is more than 6 standard deviations 
bigger than that of Ref. [M- The discrepancy increases for larger values of the cut in the 
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(MeV) 


(rcB — rLow)/rGB 


10-30 


2.0 X 10-3 


30-50 


8.7 X 10-3 


50-70 


1.9 X 10-2 


70-90 


3.3 X 10-2 


>90 


4.9 X 10-2 



Table 3: Relative differences in the rates between GB and Low. Listed are the quantities 

'-dE^^ for different bins in the photon energy. 



r^i^ fdT GB dTi,ow\ ,r. I (^^^^ 



41) V dE^ dE^ j 7 dE^ 



photon energy. 

Before attributing any significance to the predictions ([T0|), we will of course have to 
investigate the direct emission amplitude. Before doing so, we compare the rates for the 
GB amplitude (Q) with the ones for the Low amplitude (0) in the same photon energy bins 
as before. The results displayed in Table ^ show that the differences are rather small in 
all energy bins. This is due to the fact that E^^ — E^^^^ is only sensitive to the quadratic 
slope parameter c in numerically the smallest of the three parameters. Nevertheless, 
the difference between GB and Low is still bigger than the one-loop contribution to the 
direct emission amplitude to which we now turn. 



3. The full radiative amplitude is the sum of the GB amplitude (Q) and of a direct emis- 
sion amplitude E^)^: 

E^ = El^^ + E^^ . (12) 

In this paragraph we calculate the direct emission amplitude of 0{p'^). As shown in Ref. [ 
^, has the following general structure at this order: 

E^E = ^c^nterterm + J2 i^' + • (13) 

loops 

As already mentioned, there is no counterterm contribution to direct emission for i] 
7r+7r-7r°7. The (finite) loop contribution is exclusively due to diagrams of the topology 
shown in Fig. |T] where a photon should be appended to all charged lines and all vertices 
with at least two charged fields. 

The loop amplitude consists of a sum of two gauge invariant parts (for each loop 
diagram) A'^ and H^. Referring to Ref. [ ^ for details, we recall that both and 
depend only on the on-shell couplings of the vertices V^i, V2 in Fig. |I]. Those vertices have 
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X 




Figure 1: One-loop diagram for the general four-meson transition. For the radiative 
amplitude, the photon must be appended to every charged meson line and to every vertex 
with at least two charged fields. The vertices Vi, Vi are defined in Eq. (|T^. 

the general form in momentum space 

Vx = ao + aiPa-Pb + a2Pa-x 

+ as{x^ - Ml) + a^{y^ - M^) + a,{pl - Ml) + - Ml) 
V2 = bo + biPc-pd + b2Pc-x 

+ b,{x' - Ml) + b,{y' - Ml) + b,{pl - Ml) + bM - Ml) . 

The relevant on-shell coefficients for the various diagrams are collected in Table ^. We 
have included the diagrams with two neutral intermediate particles for completeness al- 
though they do not contribute to either A'^ or here. In general, is always zero in 
this case but A'^ may be non-zero depending on the assignment of particle labels]^ 

The ^ Stt couplings vanish for rriu = rrid. In contrast to Ref. [ |^, we keep 
the pion-loop contributions since they turn out to be much bigger than the kaon loops 
which we calculate in the isospin limit. The main contribution of direct emission arises in 
the interference with the GB amplitude @. The corresponding contributions to the rate, 
separately for pions and kaons, are shown in Table ^ 

The amplitude of 0{p'^) is completely dominated by the pion loops. We will explain 
the suppression of kaon loops after the discussion of vector meson exchange. Neverthe- 
less, the residual pion-loop contribution in the direct emission amplitude is quite small 
for almost all photon energies. Integrating the differential rate over the photon energy for 
i?^ > 10 MeV produces a correction to the branching ratio that is smaller than the error 
given in (|T0[) for the GB contribution only. It is even smaller than the difference between 
the rates for GB and Low amplitudes (cf. Table ^. The relative size of the loop amplitude 
increases with E^™ at the expense of decreasing rates. 

For the loop contributions to direct emission we only agree with Ref. [ ^ to the extent 

^In fact, the loop contribution to with two 7r° in the loop was missed in the calculation of Kl 
7r+7r^7r° [^. The change in the rate is numerically insignificant. 
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Table 4: On-shell coefficients of the vertices Vi, V2 defined in ( [T4l ) for the various loop 
diagrams of Fig. 1 in units of 1 /F^ and with Mf = (m^ - mu)B/{V3{M^ - Ml)). 



V{-Pb) T^aij>a) 

+ ^^^x[x)My[y) 

T^c{Pc)T^d{Pd) 


ao 




^2 


ho 


bi 


b2 


'/ — ^ VI -|- 
















-MF(3A'/2 - M2)/3 


-2Ml 





2Ml 


2 


-2 


TT TT — > TT TT 


- M^) 


U 


U 





z 


U 


— > 7r+ + 














tt'^tt^ — > H^H^ 


AM^M"^ /3 




2Mf 







—2 


T] ^ + 


4MFm2/3 




2Mf 







-2 


7] ^ + 

K-K+ TT + TT- 


M,V(2v/3) 













1 


-> 7r° + 


-M,V(2v^) 













1 


?7 — > 7r+ + 


















V 2 







-1/V2 


V2 




















/a" 
V 2 







-1/V2 


V2 



that they are small. Bramon et al. did not include the dominant pion loops and they did 
not calculate the interference with the bremsstrahlung amplitude. Taking the kaon-loop 
amplitude by itself leads of course to a tiny rate that is completely negligible [ ||] in 
comparison with the interference between the GB and the pion-loop amplitudes. 



4. Since there is no counterterm contribution to direct emission at 0{p'^) resonance ex- 
change can only enter at 0{p^). Starting from the list of 0{p^) vector and axial- vector 



couplings given in [ |T5|], we have scanned all possible contractions of the resonance fields. 
In the isospin limit (m„ = m^), the only surviving 77 7r+7r^7r°7 amplitude of this type 
is generated through the product of the following vector operators: 

Cv = hv w ( { , /r } ) + ^Ov £^.p. ( u'u^), (15) 



where we have adopted the notation of [ |T6l] for the coupling constants. For the other 
combinations of resonance terms, either the SU (2)-singlet nature of the rj field or the 
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(MeV) 


(rcB+DE - rGB)/rGB (pions) 


(rcB+DE - rGB)/rGB (kaons) 


10-30 


-1.4 X 10-4 


0.7 X 10-^ 


30-50 


2.4 X 10-4 


2.3 X 10-^ 


50-70 


2.6 X 10-3 


5.7 X 10-^ 


70-90 


9.4 X 10-3 


6.9 X 10"^ 


> 90 


3.4 X 10-2 


5.8 X 10-^ 



Table 5: Relative rate differences for the interference between GB and the one-loop con- 
tributions to direct emission. The notation is analogous to Table ^. 

resulting minimal number of pseudoscalar fields implies a vanishing contribution to the 
•T] 7r+7r-7r°7 amplitude. 

The couplings in Cy can in principle be determined from the phenomenology of 
vector meson decays. The decay rate for 7r°7 [|^] fixes = 0.037. For the second 
coupling 9v, one has to rely on models for the time being. Hidden symmetry predicts 
9v = 2hv [ [17|], which is compatible with the value Oy = 0.050 deduced from the ENJL 
model [ [1^]. We shall also assume that the field in (|T5]) describes a nonet of vector 
mesons. 

Integrating out the vector mesons in the Lagrangian ([T5[), one obtains the following 
effective Lagrangian of 0{p^) for the direct emission inr] ^ 7r+7r"7r°7: 



e _ QAihvevF'^^d^r] 



(16) 



This Lagrangian gives rise to the decay amplitude 



^DE.VMD = 3^^2^4 \P-n ■ Po9+- + Pv ' P+9-0 + Pv ■ P-9o+] (IV) 

9ij = tiP^j - tjPi 

which differs from formula (19) in [||]. 

For hv = 0.037 and 0.050 < 9v < 0.075, we find that this amplitude provides a 
contribution to direct emission that is smaller than the pion-loop amplitude, especially 
for large photon energies. Actually each of the three separate gauge-invariant terms in 
dTTD generates a contribution which is of the same order or even larger than the one from 
the pion loops. However, there is a strong destructive interference among the three terms 
which leads to the small results reported in Table |^. In the same Table we also show the 
total direct emission, obtained by summing loop and vector meson contributions. 
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(MeV) 


(rcB+DE - rGB)/rGB (VMD) 


(rcB+DE - rGB)/rGB (total) 


10-30 


0.6 X IQ-^ 


-0.7 X 10-^ 


30-50 


3.2 X 10^^ 


5.8 X 10-^ 


50-70 


7.3 X IQ-^ 


3.4 X 10-3 


70-90 


1.1 X IQ-^ 


1.0 X 10-^ 


>90 


1.0 X 10-3 


3.5 X 10-2 



Table 6: Relative rate differences for the interference between GB and direct emission: 
VMD (for hyOv = 2.8 x 10-^) and total direct emission. The notation is analogous to 
Table |. 

5. It is remarkable that the isospin conserving part of direct emission (both kaon loops 
and vector meson exchange) is smaller than the isospin violating component due to pion 
loops. In order to understand this suppression, we write the decay amplitude in the form 

= A+^{p+,P-,Po,k)g^+^ + A^o{p+,p_,po,k)g% + Ao+{p+,p_,Po,k)g^^ . (18) 

Since only two of the g^^ are linearly independent this decomposition is of course not 
unique but it will be useful in the limit niu = rud. 

The first observation is that the amplitude vanishes when the three pion momenta 
are equal. In the r] rest frame, this configuration corresponds to maximal photon energy. 
Therefore, gauge invariance alone implies that the complete amplitude and in particular 
the direct emission part is small in a region where direct emission has any chance at all 
against (generalized) bremsstrahlung. 

Let us now consider the isospin limit m„ = nid where only (part of) the direct emis- 
sion amplitude survives. In this case, isospin violation can only come from the electro- 
magnetic field which is the sum of / = and 1 spurions. Since G = —1 for 'qn^n^n'^ 
only the G = — 1, / = part of the photon contributes. Therefore, the three pions must 
be in an isosinglet combination and any two pions are in an / = 1 state. This implies that 
the amplitude (jl^) can be written in terms of a single invariant function with 

A_o{p+,P-,Po,k) = A+_{p^,po,p+,k) 
Ao+{p+,P-,Po,k) = A+_(po,P+,P-,fc) 

A+-{p-,p+,Po,k) = A+_{p+,p^,po,k) . (19) 

Both the exchange amplitude ([1^) and the kaon-loop amplitude satisfy these condi- 
tions. 
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Expressing, e.g., g% in terms of g1_, g^j^, the amplitude ( p^ can be written as 



A+-{p+,P-,Po,k) 




■+ 



■+ 



{po,P-,p+,k) g1 



iPo,P-,P+,k) g^+ 



(20) 



in the limit m„ = m^. Therefore, the amplitude is doubly suppressed for large photon 
energies: both the g^j and the two invariant amplitudes in ( pO| ) vanish in the symmetric 
configuration with p+ = P- = Po- In general, this is not the case for the explicitly isospin 
violating contributions proportional to m„ — as can be seen from Eqs. (0), @. 

The suppression of kaon-loop and vector meson exchange amplitudes is therefore 
a general feature of the amplitude in the limit rriu = ma, independently of the chiral 
order. The direct emission amplitude is strongly dominated by the pion loops. Since this 
residual pion-loop contribution is itself small compared to the dominant GB amplitude, 
the theoretical uncertainty of the total direct emission amplitude is also small and certainly 
negligible in comparison with the present experimental errors entering the GB amplitude. 

6. Our main results can be summarized as follows: 

i. The concept of generalized bremsstrahlung is very efficient in avoiding the propa- 
gation of uncertainties in the non-radiative decays to the direct emission amplitudes. 
In the case at hand, we have shown that the numerically important final state inter- 
actions in 1] ^ 7r+7r^7r° [ 0, ^ are easily incorporated in the GB amplitude. This 
allows for a very precise prediction of the radiative decay rate normalized to the 
non-radiative transition: 

B{r] ^ 7r+7r"7r°7; > 10 MeV) = (3.14±0.05) x 10"^B(r/ tt+tt^tt^) . (21) 

ii. The tttt loops dominate the direct emission amplitude even though they are isospin 
suppressed. Nevertheless, the one-loop amplitude is negligible compared to GB for 
most of the photon energy range. 

iii. Isospin conserving contributions to direct emission such as kaon loops or vector 
meson exchange are suppressed, especially for large photon energies. Since this is 
a general feature to all orders in the chiral expansion we expect it to be very difficult 
if not impossible to observe any direct emission effect in r/ ^ 7r+7r~7r°7 even with 
the anticipated yield of 10^ i] per year [ pT|]. 
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